Influence of Noise on Force Measurements 
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We demonstrate how the ineluctable presence of thermal noise alters the measurement of forces 
acting on microscopic and nanoscopic objects. We quantify this effect exemplarily for a Brownian 
particle near a wall subjected to gravitational and electrostatic forces. Our results demonstrate that 
the force measurement process is prone to artifacts if the noise is not correctly taken into account. 



The concept of force plays a central role in our un- 
derstanding of nature. Due to the ongoing trend to- 
wards miniaturization, the investigation of forces relevant 
at microscopic and nanoscopic length scales is attract- 
ing an increasing amount of attention. Examples range 
from the elastic properties of biomolecules [jy] to Casimir 
forces Instrumental in such trend has been the in- 
vention of new methods to measure ultra-small forces in 
a range down to few femtonewtons 0, 0]. Apart from 
the technological challenge intrinsic to measuring such 
minute forces, it is important to realize that the general 
concept of how forces are measured in macroscopic sys- 
tems cannot be simply scaled down to nanoscopic objects, 
mainly due to the presence of thermal noise affecting the 
motion of small objects (Brownian motion). As will be 
demonstrated below, force measurements in the presence 
of thermal noise are prone to artifacts, unless the noise is 
properly taken into account. However, despite the great 
number of experiments measuring forces in microscopic 
systems, the role of noise has not been adequately ad- 
dressed yet. 

In this Letter, we measure the forces acting on a Brow- 
nian particle in front of a wall: using two widely em- 
ployed force-measurement methods, we obtain strongly 
contrasting forces, which deviate in their magnitude and 
even their sign. We track this disagreement down to a 
well-known mathematical controversy about the interpre- 
tation of stochastic differential equations in the presence 
of a diffusion gradient and we demonstrate how exper- 
imental data should be analyzed in order to obtain the 
correct underlying forces and to avoid artefacts. 

When a microscopic body is suspended in a liquid, vis- 
cous forces prevail by several orders of magnitude over in- 
ertial effects Q . The presence of a constant external force 
F in such an overdamped system results in a constant 
particle drift velocity — F/^, where 7 is the particle's 
friction coefficient. Since Vd — dz/dt can be retrieved 
from the measured particle displacement dz within time 
dt, the force can accordingly be measured as F = ^Vd- In 
the case of large forces this obviously leads to a univocal 
result. However, when the force exerted on the object is 
comparable or even smaller than the random forces due 
to the permanent collisions with the surrounding liquid 
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FIG. 1: (a) A Brownian particle (drawn not to scale) diffuses 
near a wall in the presence of gravitional and electrostatic 
forces. Its trajectory perpendicular to the wall is measured 
with TIRM. (b) Comparison of measured (bullets) and calcu- 
lated (line) vertical diffusion coefficient as a function of the 
particle-wall distance, (c) Experimentally determined prob- 
ability distribution of the local drift dz for dt — 5 ms at 
z = 380 nm (grey). The dashed line is a gaussian in excellent 
agreement with the experimental data. 



molecules, a different value of the force is measured in 
each experiment. Typically one can average over many 
independent measurements 



F = jVd, 



(1) 



{dzj) and j denotes the j-ih 



where we define Vd = 
experimental value. 

Indeed Eq. ([T]) is key in non-equilibrium force mea- 
surements, e.g. studying the relaxation of biomolecules 
0-11] , colloidal interactions 0, ITol. or Brownian motion 
in non-conservative force fields jllT - flsj . It must be em- 
phasized, however, that Eq. ([1]) is only valid for a spa- 
tially homogeneous diffusion coefficient D of the object to 
which the force is applied; within linear response this is 
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FIG. 2: Comparison of forces obtained from local drift veloc- 
ity measurements of a Brownian particle according to Eq. ([T]) 
(bullets) and measured from the equilibrium particle height 
distribution (squares). The solid line corresponds to a fit to 
Eq. (HJ. 



related to its friction hy D — UbT/^ 13]. Generally, Eq. 
(HI) does not hold when D varies vi^itli position, e.g. due 
to hydrodynamic interactions betvifeen the particle and 
nearby walls or other particles, a situation often encoun- 
tered in experiments. It has been shovi^n that such spatial 
variations in the diffusion coefficient have to be explicitly 
considered when numerically computing the particle tra- 
jectory IJ, |l5| . As a consequence, application of Eq. ^ 
leads to erroneous forces, which may severely affect the 
physical interpretations of experimental data. 

To demonstrate the effect of noise on force measure- 
ments, we experimentally study a colloidal particle (ra- 
dius R = 1.31 ± 0.04 /im, density pp = 1.51 g/cm^, MF- 
F-1.3, Microparticles GmbH) immersed in water (density 
Ps = 1.00 g/cm^) and diffusing in a closed sample cell 
above a planar wall, placed at z = (Fig. [Ija)). This is 
arguably the simplest realization of a diffusion gradient. 

The particle's trajectory perpendicular to the wall z[t) 
(Fig. [TJa)) is sampled with nanometer resolution at a 
sampling rate of 1 kHz over 200 minutes employing a sin- 



gle particle evanescent light scattering technique known 
as total internal reflection microscopy (TIRM) [ll,[l3l. A 
p-polarized laser beam (A = 658 nm) is totally internally 
reflected at a glass-liquid interface generating an evanes- 
cent field decaying into the liquid. The particles trajec- 
tory is obtained from the scattering intensities, which 
depends on its position relative to the interface. 

All conservative forces acting on the particle are [3l 



F{z) = Be- 



G. 



eff- 



(2) 



The first term is due to double layer particle- wall forces, 
with = 18 nm the Debye length (300 p,M NaCl salt) 
and B a prefactor depending on the surface charge den- 
sities. The second term accounts for the effective gravi- 
tational contribution Ge// = ^TTR^{pp — Ps)9, with g the 
gravitational acceleration constant. 

Far away from any surface, the diffusion coefficient 
of a spherical particle is Doo = ksT /QiTrjR^ where 77 
is the liquid shear viscosity. Close to a wall, however, 
the diffusion coefficient sharply decreases due to hydro- 
dynamic interactions. From the solution of the corre- 
sponding creeping flow equations one obtains an analyt- 
ical expression for D \ (z), the component of D perpen- 
dicular to the wall [18|, which is plotted for our exper- 
imental conditions as solid line in Fig. [TJb). The cor- 
responding data obtained from the experimentally mea- 
sured particle trajectory according to the conditional av- 
erage Da^{z) ^^^{[z{t + dt)-z{t)Y I z{t) = z) [ll,[20] 
(symbols in Fig. (Tt) show excellent agreement with the 
theoretical prediction; due to the particle-wall electro- 
static repulsion only distances above 180 nm are sampled. 

Since in our system the force depends on z, the average 
drift velocity Vd in Eq. ([T|) has to be replaced by its local 
value Vdfyz') = {z{t -\- dt) — z{t) \ z{t) = z). The time 
interval dt for which the displacement is considered has 
to be sufficiently small to guarantee that the force act- 
ing on the particle can be treated as locally constant. In 
our experiments this condition is met for dt < 10ms. Fig. 
[ijc) shows the probability distribution of Vd for dt = h ms 
and z = 380 nm, which almost perfectly agrees with a 
gaussian distribution and thus confirms that within such 
small time steps, the spatial variation of the force can be 
neglected. After having replaced the constant friction co- 
efficient in Eq. ([T]) with its local value 7(2;) = ksT / D±{z) 
one finally obtains the local force F{z) acting on the par- 
ticle. The result is shown as bullets in Fig. [S] 

Since our system is in thermal equilibrium, the forces 
acting on the particle can be also obtained from the 
measured particle-wall interaction potential U (z) = 
—kBT\np{z)^ where p{z) is the experimental particle's 
position equilibrium distribution. Contrary to F{z), this 
approach is valid independently of the additional pres- 
ence of hydrodynamic interactions. The corresponding 
force F[z) = —-^U{z) is shown as squares in Fig. [5] 
and systematic deviations from F{z) are evident. In ad- 
dition, the force-distance relation obtained via the p{z) 
is in quantitative agreement with Eq. ([2]) (solid line in 
Fig. [2]), where Ge// and are taken from the exper- 
imentally known parameters. The prefactor B has been 
treated as a fit parameter and its value B — 770 pN is 
in good agreement with other TIRM experiments under 
similar conditions 16|. In the following, we discuss the 
reason for this discrepancy providing a simple method 
to correctly interpret the results of non-equilibrium mea- 
surements and to reliably measure forces. 

The motion of a Brownian particle can be described 
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by a stochastic differential equation (SDE) where a ran- 
dom function is added to an ordinary differential equa- 
tion (ODE) [2lj. This approach was introduced at the 
beginning of the 20"^ century by Smoluchowski, Einstein, 
Langevin, and Kolmogorov 22-241 and put on a firmer 
mathematical ground in the 1950s and 1960s by Ito and 



Stratonovich |25ll26l |. For a Brownian particle in the pres- 
ence of a variable diffusion coefficient D±{z), e.g. near 
to a wall, the corresponding SDE reads 



dz 



Fjz) 
7(z) 



dt + ^2D^{z)dW, 



(3) 



where is a Wiener process, i.e. a stochastic process 
almost surely continuous, almost nowhere differentiable, 
and whose increments dW are stationary, independent 

Integration of Eq. ([3]) 



and normally distributed 
yields 



211 



r ^dt+ r ^/mii^dw. (4) 

Jo 7^ Jo 



z{T) = z(0) + 



Due to the irregularity of the Wiener process the value 
of the stochastic integral on the r.h.s. is ambiguous. It 
is defined as the limit of integral sums where the inte- 
grand is evaluated inside each bin at a given position 
parametrized by a € [0,1], i.e. jj" yj2D^{z)dW — 

lim^^oo Y.n=o V2^J-(^(^n))AW„, where t„ = ^^T. 
Since 14^ is a function of unbounded variation, differently 
from ordinary Riemann-Stieltjes integrals, this leads to 
different values for each choice of a. A loose understand- 
ing of such indetermination can be gained by considering 
14^ as a random sequence of pulses, each having an in- 
finitesimal duration, but still a finite amplitude; the value 
of a determines at which position during each jump the 
integrand should be evaluated [i^l • Common choices are 
a — Q (the Ito integral)^ a = 0.5 (Stratonovitch integral), 
and a = I {anti-Ito or isothermal integral). This is in 
sharp contrast to ODEs, which have a univocal interpre- 
tation. 

The values of stochastic integrals for different a are 
related to each other by a precise mathematical relation- 
ship [2§| . For example, the solution of Eq. dU can always 
be written as an Ito integral to which an a-dependent 
correction term is added 
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FIG. 3: Calculated noise-induced drift (lines) for different 
values of a. The symbols correspond to the experimentally 
determined noise-induced drift. 



the latter only disappears when the diffusion coefficient 
is homogeneous; otherwise, it has to be accounted for in 
deducing the forces ^"(2;) acting on a particle from its 
measured drift velocity Vd{z). We remark that, while 
an equilibrium measurement of the potential constitutes 
a unambiguous means to experimentally determine the 
forces, in non-equilibrium situations it is not a priori 
clear which value of a should be used and therefore it 
is not a priori clear how such noise-induced drift has to 
be considered in non-equilibrium experiments. 

From the theoretical D^{z) 0] one immediately ob- 
tains the noise- induced drift a '^"^^J^^ , which is plotted 
as lines for a = 0, 0.5, 1 in Fig. [3] In order to deter- 
mine which value of a is valid in our experiment, we 
have calculated the noise-induced drift from the experi- 
mentally measured forces. According to Eq. ^ it is de- 
termined from the difference of the forces obtained from 
the drift-velocity measurement and from the equilibrium 
potential measurement, both presented in Fig. [21 i.e. 
{F{z) — F{z))/j{z). The experimental data (symbols 



F(z) 



z{T) = z(0)-f / ^^dt+ / ^/2DAz)dW 
lo l[z) Jo 



+a 



dD^{zf 
— ail 



dz 



From this one obtains the particle drift velocity 

"f{z) 7(z) dz 



(5) 



(6) 



in Fig. show good agreement with the noise-induced 
itift obtained for a — 1, i.e. the isothermal integral. 
All other choices of a, in particular negligence of noise- 
induced drift (a = 0), lead to significant differences. 
In Fig. |4]we have plotted as open symbols 



F{z) — q;7(z) 



dD^{z) 
dz ' 



(7) 



The first term on the r.h.s. is the deterministic drift 
due to the "real" forces acting on the particle, while the 
second term represents a noise-induced drift. Obviously, 



i.e. the experimental F{z) (open symbols in Fig. [2]) 
with a correction term for a = 1 (squares), a = 0.5 
(triangles), and a = (bullets), respectively. It should 
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FIG. 4: Forces obtained from a drift velocity experiment with 
added noise-induced drift (see Eq. ([7]) with a = 1 (open 
squares), a = 0.5 (open triangles), and a = (open dots). 
The solid squares represent the forces obtained from an equi- 
librium measurement (same as in Fig. [2}. 



be emphasized that not only the absolute value but even 
the sign of the force depends on the choice of a. As 
closed symbols we have superimposed the forces F{z) 
obtained from the particle probability distribution (same 
as solid symbols in Fig. [2|) , which shows good agreement 
for a = 1. Since the gradient of the diffusion coefficient 
vanishes far away from the surface, the force-dependence 
on a is most pronounced close to the wall but weakens 
at larger z. 

Our experiments clearly demonstrate that the pres- 
ence of noise-induced drift has to be considered in non- 
equilibrium force measurements; otherwise, this can lead 
to artifacts in the measured forces, which may even sug- 
gest the wrong sign of the force. While the correction is in 
our case on the order of several femtonewton, it becomes 
more significant in the presence of larger diffusion gradi- 
ents, i.e. for shorter particle- wall distances or for smaller 
particles. We stress, furthermore, that a constant dif- 
fusion can be assumed only for a particle far from any 
boundary. Such boundaries are naturally introduced by 
surfaces or by other particles in suspension, a situation 
that is typically met in experiments. 

In the cases in which thermodynamic consistency must 
be satisfied, a = 1 is the correct choice [2^, [s^l- This is 
particularly true for the experimental system we have in- 
vestigated, i.e. a colloidal particle performing Brownian 
motion coupled to a thermal bath. Indeed, the conven- 
tion Of = 1 is the only one naturally leading to the usually 
accepted steady-state probability distribution; other con- 



ventions require the addition of a spurious drift term to 
the L ang evin equation to account for the noise-induced 
drift [30] ■ One is also naturally led to the convention 
a — 1 when considering the vanishing mass limit of the 
second-order Smolukowski equation, which, despite con- 
taining a random term, has an unambiguous interpre- 
tation even in the presence of a diffusion dependent on 
position, but not on velocity. We remark, however, that 
in case of a velocity-dependent diffusion coefficient also 
the interpretation of the Smolukowski equation becomes 
ambiguous |31i] . 

In a more general sense the value of a depends on the 
system under study [s^]- What works for the motion 
of a Brownian particle might not be appropriate for the 
description of other stochastic processes, e.g. stock mar- 
ket behavior or ecosystem dynamics. This study does 
not claim to find out which is the "right" value of a 
for all the situations modeled by SDEs. Nonetheless, it 
demonstrates that the intrinsic ambiguity of SDEs with 
multiplicative noise is indeed amenable to experimental 
scrutiny and it may clear the way for similar studies in 
other fields that make an intensive use of SDEs, such as 
economics and biology. 
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